submitted to ApJ 



Suzaku Observations of the Circinus galaxy 

Y. Yang 1 ' 2 , A. S. Wilson 1 , G. Matt 3 , Y. Terashima 4 , L. J. Greenhill 5 

ABSTRACT 

We report Suzaku observations of the active, Compton-thick Circinus galaxy. 
Observations were obtained with both the X-ray Imaging spectrometer (XIS) 
and the Hard X-ray Detector (HXD). Below 10 keV, the nuclear spectrum is 
dominated by radiation reflected from cold dense gas of high column density, 
while above 13 keV the radiation is directly transmitted nuclear emission seen 
through a column density of ~ 4 x 10 24 cm~ 2 . In the 0.2-10 keV band, the 
XIS spectrum is contaminated at 5% level by the brightest off-nuclear source 
in Circinus (CG X-l), but drops to 1% in the 5-10 keV and is negligible at 
higher energies. We find no significant evidence for variability in the hard (> 
12 keV) emission. The Circinus is marginally detected with the HXD/GSO in 
the 50-100 keV band at 2.5a level. We model the 3-70 keV band XIS+PIN 
spectra with a four components: the Compton transmitted nuclear emission, 
the reflected nuclear emission, a soft power law (representing a combination of 
scattered nuclear emission, extended emission and contamination by sources in 
the galaxy below a few keV). The hard nuclear power- law is found to have a 
photon index I\ ~ 1.6, very similar to the soft power-law. The high energy 
cut-off is Ec — 49 keV. These results agree with those from BeppoSax. An 
extrapolation of this model up to the GSO band shows good agreement with the 
GSO spectrum and supports our detection of the Circinus up to ~ 100 keV. 

Subject headings: galaxies: individual: Circinus Galaxy — galaxies: Seyfert — galaxies: active- 
X-rays:galaxies 
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1. Introduction 

The Circinus galaxy is a large, nearby (4 ± 1 Mpc; Freeman et al. 1977) galaxy that 
harbors both a circumnuclear star burst and a Seyfert 2 nucleus. Evidence for an obscured 
Seyfert 1 nucleus is provided by the finding of a broad (FWHM ~ 3000 kms" 1 ) Ha line 
component in polarized light (Oliva et al. 1998). The picture of an obscured Seyfert 1 nucleus 
is also supported by both the discovery of highly ionized gas extending along the minor axis 
of the galaxy, with a morphology that is reminiscent of the ionization cones seen in other 
Seyfert galaxies (Marconi et al. 1994; Wilson et al. 2000), and direct X-ray detection of the 
nucleus through a column density of ~ 4 x 10 24 cm -2 (Matt et al. 1999, hereafter M99; Soldi 
et al. 2005). 

ASCA, BeppoSax and XMM-Newton spectra below 10 keV show both a flat (hard) 
spectrum and high equivalent width Fe Ka emission, characteristic of Compton reflection 
from cold gas illuminated by a power-law continuum. A single scattering Compton shoulder 
was found in Chandra HETG observations (Bianchi et al. 2002) and was confirmed by 
XMM-Newton (Molendi, Bianchi, & Matt 2003), demonstrating conclusively the existence of 
Compton thick matter. Because of their very flat hard X-ray spectra, Compton-thick AGNs 
are a key ingredient in the synthetic model of the hard cosmic X-ray background (Setti & 
Woltjer 1989; Madau, Ghisellini, & Fabian 1994), which peaks at 30-40 keV and rolls over 
at higher energies. After NGC 4945 and NGC 1068, Circinus galaxy is the third brightest 
Compton-thick AGN that allows detailed study. 

In this paper, we present a broad band observation of the Circinus galaxy using the 
Suzaku telescope. The low background and high sensitivity of Suzaku allow improved con- 
straints on the nuclear continuum and the nature of the absorbing and reflecting matter. In 
§ El we describe the observation and data reduction. The results are presented in § [31 We 
discuss the results and their implications in § |U 

2. The Suzaku Observation and Data Analysis 

The Circinus galaxy was observed for ~ 100 ks with Suzaku on July 21, 2006. The data 
were obtained with the four X-ray Imaging Spectrometer (XIS,Koyama et al. 2007) CCDs, 
and the non-imaging, collimated hard X-ray detector (HXD). Each of the four XIS CCDs 
(XIS0, 1, 2 and 3) is located in the focal plane of a foil X-ray telescope (Serlemitsos et al. 
2007). The CCD chip for each camera has a dimension of 1024 x 1024 picture elements 
("pixels"), and covers a 17.8' x 17.8' region on the sky. The CCD cameras are sensitive in 
the range 0.2-12 keV. One of the XISs, XIS1, uses a back-side illuminated (BI) CCD, while 
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the other three use front-side illuminated (FI) CCDs. The BI CCD has better quantum 
detection efficiency at sub keV energies. Each CCD sensor has two 55 Fe calibration sources 
located on the side wall of the housing; these sources illuminate two adjacent corners of the 
CCD. The spectral resolution of the XIS at 6 keV is 130 keV. The half power diameter of 
the X-ray telescope is 2'. The HXD uses a novel well-type phoswich counter (Takahashi et 
al. 2006) which greatly reduces the instrument background. The HXD sensor consists of 
16 main detectors and surrounding 20 crystal scintillators for active shielding. Each unit 
consists of two types of detectors: a GSO scintillation counter, and 2 mm-thick PIN silicon 
diodes located inside the well and in front of the GSO scintillator. The PIN diodes are most 
sensitive below ~ 60 keV, while the GSO is sensitive above ~ 40 keV. Our observations were 
performed during the period when bias voltages for 16 out of 64 PIN diodes were reduced 
from 500 V to 400 V to suppress the rapid increase of noise events that possibly result from 
in-orbit radiation damage. Since the thickness of the depletion layer depends on the bias, 
this change slightly affected the energy response of the 16 PIN diodes especially in the energy 
band higher than 20 keV. 

The observation started on July 21, 2006 at 12:29:57 UT, and stopped on July 23, 2006 
at 1:08:24 UT. The nominal point of the observation lying at the center of the XIS field, is 
at «2ooo — 14 /l 13 m 09 s .84, S 2 ooo — — 63 o 20'24". The data were first reprocessed with version 
1.2.2.3 calibrations and analyzed with FTOOLS included in HEADAS 6.2. Preliminary 
results were presented in (Yang et al. 2007). After we finished the first draft of the paper, 
version 2 reprocessed data became available. Improvements of the XIS calibration have 
resolved some of the issues in the XIS spectra we found when using version 1 processing. We 
therefore use the version 2 reprocessed XIS data (process version 2.0.6.13). For the HXD 
data, we found no significant difference in our PIN spectra using either version of data. 
Since the reproductivity of version 2 processed HXD/GSO background is still unavailable, 
which affects our GSO results, we opt to keep using version 1 processed HXD data, where 
the systematics are better known. Following the recommendation of the Suzaku instrument 
team, we selected only the 48 PIN diodes that used a bias voltage of 500 V. The cleaned 
exposure times for XIS and HXD are 108.3 ks and 86.0 ks, respectively. 

We extracted the XIS spectra and light curves from a 2.5' radius region centered on 
the Circinus galaxy. The background spectra were extracted from a set of circular regions 
away from the Circinus galaxy and other bright point sources and extended emissions. The 
spectra from the calibration sources were also extracted. The broad point spread function of 
the XRTs mixes the emission from the nucleus with the emission from several point sources, 
including the two bright ultra-luminous X-ray sources (ULXs), named CG X-l, CG X-2 in 
Bauer et al. (2001), and that from the diffuse gas seen in Chandra images (Smith & Wilson 
2001). This may complicate the interpretation of the soft X-ray spectrum. We show the 
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XIS extract region overlaid on the Chandra image of the Circinus galaxy field in Fig. [T] to 
illustrate the location of the ULXs. 

Below 100 keV, the passive collimator of the HXD defines a 34' x 34' field of view. 
Instead of rocking between source and blank sky, the non-cosmic background of the HXD 
is modeled based on the earth occultation data, monitoring data of the actual particle flux, 
and the orbit data of the satellite. The accuracy of the background is determined by the 
reproductivity of the model rather than the statistical error of the background counts. Since 
the non X-ray background contributes significantly to the HXD counts, the accuracy of the 
model background is the major factor determining the sensitivity of the instrument. The 
accuracy of the model background can, in most cases, be tested by comparing the simulated 
background with that obtained during Earth occultation. This, however, is not the case 
for the Circinus galaxy because the galaxy is located in the continuous viewing zone of 
Suzaku with earth elevation angle always > —5°. We therefore rely solely on the HXD non 
X-ray background files provided by the HXD team in our analysis. The uncertainty of the 
background can only be stated in a statistical sense. 

3. Results 
3.1. Variability 

3.1.1. XIS 

The brightest (L x = 3.71 x 10 39 ergs -1 in the 0.5 to 10 keV band, corrected for absorp- 
tion) off nuclear source in Circinus - CG X-l - was found to be periodic with period 27.0 ks 
through Chandra observations by Bauer et al. (2001). We examine the folded 0.2-10 keV 
band XIS light curve (Fig. [2]) and the period is clearly detected. The variation amplitude 
between minimum and maximum is 10%. However, the amplitude reduces to only 2% in the 
5-10 keV band, indicating the spectrum of the variable source should be steeper than that 
of the nuclear source. This result agrees with the notion that the spectrum in this band is 
dominated by the reflected component of the nuclear emission. 

3.1.2. PIN 

On the other hand, the observed 12-50 keV PIN light curve suggests a slow variability 
on a time scale of ~ 50 ks during our observation, even after the time variable deadtime 
correction has been performed. The standard deviation of the count rate is 15% (4 ks 
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binning). However, the model background count rate is ~ 60% of the total count rate and 
thus can be an important source of the observed variability. For 4 ks bins, the uncertainty 
of the PIN non X-ray background reproductivity measured by the standard deviation of the 
residual between the earth occultation and blank sky data and the background model is 
o ~ 7% (Mizuno et al. 2006, Fig. 3). The residual also tends to be larger during the SAA 
passages (because the background is harder to model during these periods), and produces 
a slow variation on time scales < 70 ks. This level of variability translates to a a ~ 10% 
uncertainty in the source count rate, which is about the same level of variation in our 
PIN light curve. Therefore, our detected variability in the hard X-ray emission is probably 
not real. This conclusion is consistent with the emission being Compton scattered in a 
thick torus and the intrinsic variability has been suppressed and smoothed by the scattering 
medium. The lack of variability in Circinus galaxy compares interestingly to NGC 4945, 
which show strong variability in hard X-ray (Iwasawa et al, 1993; Madjeski et al. 2000 and 
Itoh et al. 2008). NGC 4945 has very similar spectral properties as Circinus galaxy with 
N H ~ 5 x 10 24 cm 2 (Itoh et al. 2008). Using Monte Carlo simulations, Madjeski et al. 
(2000) show that the large variability is consistent with the absorber in NGC 4945 being a 
geometrically thin disk. Since Circinus galaxy has similar optical depth as NGC 4945, we 
can estimate the geometry of the torus by applying the same simulation result shown in Fig. 3 
of Madjeski et al. (2000). Assuming a Fe abundance of unity, our upper limit of variability 
(< 15%) suggests the half angle subtended by the torus viewed at the central source > 80°. 
This is roughly consistent with the estimate subtended of the opening angle of the torus in 
Circinus galaxy by Ghisellini et al. (1994). 



3.2. GSO detection 

The background subtracted GSO count rate for the Circinus galaxy in the 50-100 keV 
band is 0.11 cts/s. The corresponding count rate of the background is 5.43 cts/s. The 
uncertainty of the background model for a 1 day exposure in the 50-100 keV band is ~ 0.81% 
(Takahashi et al. 2007, Fig. 3). Therefore, the Circinus galaxy is detected in this band at a 
~ 2.5(7 level. In Fig. [3] we show the extrapolation of the 3-70 keV XIS+PIN spectra to the 
GSO band up to 200 keV using the spectral model described in § 13.3.31 A good agreement is 
found between the data and the extrapolation of the model. Better constraints on the GSO 
spectrum await the improved accuracy of the GSO background. 
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3.3. The Spectra 

In Fig. [3] we show the Suzaku spectrum in the 3-200 keV band. Before modeling the 
broad band spectrum of Circinus Galaxy, it is instructive to first analyze the XIS and HXD 
spectra separately. This allows us to understand the calibration of the instruments and the 
contamination from off nuclear sources in each instrument. Below 10 keV, the spectrum is 
dominated by reflected nuclear emission, while above 20 keV, the spectrum is almost com- 
pletely Compton transmitted emission. Analyzing the XIS and HXD spectra separately make 
it easier to understand each component of the spectrum without too much confusion from 
other components. Degeneracies between parameters such as the intrinsic photon index and 
the absorption column density are inevitable when using only the XIS or PIN spectrum, and 
are best addressed with joint analysis using both XIS and HXD spectra, which is discussed 
in 5 13X31 

3.3.1. The X-ray Imaging Spectrometer Spectra 

Within the 2.5' extraction radius, the most likely contamination of the XIS spectrum 
comes from the two ULXs in the field (Bauer et al. 2001; Smith & Wilson 2001). Without 
simultaneous Chandra observations, it is hard to estimate the true contribution from these 
off-nuclear sources. In most of the previous observations, the 2-10 keV fluxes of the two 
ULXs remain about an order of magnitude lower than that of the nucleus. However, on 
March 14, 2001, the Chandra flux from CG X-l reached the same level of the nuclear flux 
(Smith & Wilson 2001; Bianchi et al. 2002). We assume that the contamination from the 
ULXs and other unresolved soft emissions can be represented with a single power-law. This 
component dominates the soft X-ray emission below a few keV, and we call it the soft power- 
law component in the following discussions. By modeling the XIS spectrum with a reflection 
component from dense cold gas near the blackhole (pexrav model in XSPEC, Magdziarz 
& Zdziarski 1995) plus the soft power-law with ad hoc added emission lines, we found the 
total 2-10 keV flux within the extraction region to be 1.73 x 10~ n ergcm~ 2 s -1 . This flux 
lies between that of the two BeppoSax observations on Mar. 13, 1998 and Jan. 7, 2001. 
The soft power-law component is found to have T s = 1.54 and Nh = 3.8 x 10 21 cm -2 , 
where the absorption column density includes both that from the Circinus galaxy and the 
our own Galaxy. The 2-10 keV flux of the component is 5.76 x 10 -12 erg cm -2 s" 1 . This is 
consistent with the best-fit soft power- law in M99, but higher than the average flux of either 
CG X-l (9 x 10~ 13 erg cm" 2 s" 1 ) or CG X-2 (1.3 x 10" 12 erg cm" 2 s" 1 ). This suggests that 
1) the ULXs are likely to be in states similar to those during the March 1998 observation, 
and 2) the soft emission may come mostly from the scattered nuclear emission rather than 
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from the ULXs. Moreover, since the photon index of the soft power-law is higher than 
that of the cold reflection component in our model, the relative contamination from the 
ULXs should drop significantly at higher energies. In the 5-10 keV band, the soft power-law 
component contributes only ~ 20% to the total flux. The very low amplitude of the 27 ks 
periodical variation at high energies also suggests that the spectrum of CG X-l is steep 
in this observation. Nevertheless, the contribution from the ULXs may still be significant 
sources of contamination in the energies below 2 keV. For this reason, we focus our analysis 
on the hard X-ray spectrum > 3 keV, and defer the detailed study of the soft X-ray emission 
to a future paper. 

We compare the XIS energy scales by fitting the strongest emission lines in the 3- 
10 keV spectra. Small, but significant, differences are seen in the line energies between 
different XIS spectra. The largest difference is found between the line energies using the BI 
XIS-1 spectrum and those using the FI CCD spectra. The energies of lines in the XIS-1 
spectrum are higher than those from the FI CCDs. For example, the Fe Ka line in the 
XIS-1 spectrum is ~ 25 eV higher than the best-fit energy of the line using the rest of the 
XIS detectors. This difference is significantly larger than the nominal 0.2% uncertainty level 
of the energy calibration. On the other hand, such discrepancy is not seen in the emission 
lines from the calibration sources located at the corners of the CCDs. This problem has 
been seen in observations taken after mid-2006 in observations without spaced-row charge 
injection (Hamaguchi 2007, private communication). A possible cause is that the charge 
transfer inefficiency is not properly corrected in these observations. The best-fit Fe Ka line 
rest frame energy using the FI CCDs is 6.401 keV, which is very close to the theoretical 
value of the doublet (6.391 keV and 6.405 keV). This is also the case for the rest of the 
emission lines in the 3-10 keV spectrum, indicating the energy calibration of the FI CCDs 
in our observation are probably better than that of the BI CCD. Given that the effective 
area of XIS-1 is significantly lower in the 3-10 keV band than the FI CCDs, dropping the 
data from XIS-1 does not significantly change the signal-to-noise, but significantly improves 
the statistics of the fit. Thus, we only analyze the spectra from the FI CCDs in this paper. 

The model we employed to fit the XIS spectra consists of three components: the reflected 
nuclear emission (pexrav in XSPEC), a soft power-law that represents the scattered nuclear 
emission and the emission from contamination sources, and a set of emission lines. To 
reduce the number of free parameters, we have adopted in the pexrav model the inclination 
angle of cos(z) = 0.45 and set the reflection parameter to -1 to make the component pure 
reflection. Since the normalization and the inclination of the slab i is highly degenerate, 
choosing cos(i) = 0.45 (or i ~ 63°) is mainly to reduce the number of free parameters. For 
a cylindrical torus, the reflecting surface is the inner side of torus which is perpendicular 
to the accretion disc. Therefore, choosing % = 63° means we are looking at the torus at 
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27 °. This choice is consistent with the estimates that the inclination of the system is < 40° 
(Matt et al. 1999). We also assume the reflecting gas to have the same metallicity as that 
of the Sun except for element Fe, which is left as a free parameter^. The upper cutoff 
energy of the spectrum of the reflected component is poorly constrained with only the XIS 
spectrum, and the value from the PIN spectrum Ec = 48.7 keV (see § 13.3.21) is adopted. 
We consider the following two cases for the continuum: 1) The soft power-law parameters 
are frozen to those found in the 0.3-10 keV spectral fitting. This component can be a 
combination of emissions from off-nuclear sources and the scattered emission. The intrinsic 
photon index of the reflected component is found to be 1.791^3; 2) If the emission above 3 
keV is predominantly the nuclear emission scattered by ionized matter thus the photon index 
of the soft power-law and the intrinsic photon index of the reflection component should be 
the same (with the assumption that scattering dominates over the line emission, which is 
justified above 3 keV), and are thus linked in the spectral fitting. The photon index in this 
case is found to be 1.56jV,?. 



3.3.2. The Hard X-ray Detector spectrum 

Besides the two bright ULXs, CG X-l and X-2, another bright point source located 
at ~ 5' SW of the Circinus galaxy, is also within the field of view of the HXD and could 
be a source of contamination. The source can be seen in the full band XIS image, but 
not detected in the 10-13 keV image, suggesting the contribution from this source is weak 
compared to the nucleus of Circinus galaxy in the hard band and thus can be ignored in the 
HXD spectrum. 

We do not include the PIN data below 13 keV because of the considerable noise in 
this band, likely caused by thermal events. Since the Circinus galaxy is only ~ 2.5a above 
the GSO background in the 50-100 keV band, the uncertainty in the GSO background can 
introduce significant uncertainty in the spectrum. We therefore choose not to include the 
GSO data in the spectral modeling. This, however, limits our ability to constrain the high 
energy cut-off of the nuclear emission, which is likely to be ~ 50 keV. 

The cosmic background in the PIN band is estimated using the HEAO-1 best-fit spec- 
trum (Boldt 1987). 

(£j \ -0.29 / E \ 
-r^7 exp — - erg cm -2 s" 1 str^ 1 keV" 1 
3 keV / \ 40 keV / 



6 We use the metallicity table from Anders & Grevesse (1989) 
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and the response file for uniform emission in an area of 2° x 2°. We use XSPEC to simulate 
the CXB spectrum and subtract it from the PIN spectrum. The contribution of the CXB 
to the PIN count rate is ~ 5.7%. The Circinus galaxy lies ~ 10° away from the the bulk of 
the Galacic ridge emission (Warwick et al. 1985). Using the spatial and spectral models of 
the Galactic ridge emission (Yamasaki et al. 1997; Krivonos et al. 2007), we estimated the 
PIN count rate from the Galactic ridge emission to be < 1% and thus ignored. 

The PIN spectrum above 20 keV is well described by a heavily absorbed power-law with 
high energy cut-off, which probably originates from the nuclear emission transmitted through 
a torus or thick disk. We fit the 20-70 keV PIN spectrum with a Compton transmission 
model based on Monte Carlo simulations in a spherical distribution of matter (see M99 and 
references therein). The model includes Compton down scattering and the decline of the 
Klein-Nishina cross section at high energies. A solution with photon index T = 1.89lo32^o 11 
and E c = 48.71 

^-il keV is found. The errors are the statistical errors that correspond to 
the 90% confidence ranges (the first numbers), and the systematic errors caused by the ~ 4% 
(90% confidence) uncertainty of the PIN background (the second numbers). We estimate 
the systematic errors by varying the normalization of the PIN background by ±4%. In the 
following, if only one set of errors are quoted, they correspond to only the statistical errors. 
The photon index and the high energy cut-off values agree very well with the BeppoSax 
results. The photon index also agrees with the best-fit XIS value. However, solutions with 
flatter spectra and lower cutoff energies also exist. In Fig. HI we plot a contours map of x 2 as a 
function of T and Eq- The plot shows multiple local minima with overlapping 2a contours. 
In fact, some models with very flat spectra (r < 1) produce slightly but insignificantly 
smaller x 2 values. This makes XSPEC converge to unphysical solutions when the parameter 
space is searched. We therefore fix Ec = 48.7 keV when studying the confidence range 
of other spectral parameters (§ 13.3.31) . The best-fit column density is found to be Nh = 
5.OI2 x 10 24 cm -2 , which agrees very well with previous results. The uncertainty of Nh 
is large because the parameter is best constrained in 10-20 keV band. However, the reflection 
component also contributes significantly in this band and thus Nh is best determined by joint 
fit of the XIS and HXD spectra. 



3. 3. 3. Overall 3-70 ke V spectrum 

We now study the broad band spectrum using both the XIS and HXD spectra. Again 
only the XIS spectra from 3-10 keV and the HXD/PIN spectrum between 13-70 keV were 
used in our analysis. Following M99, the broad band emission is modeled by combining the 
XIS and PIN models we have described. In this model, the nucleus is obscured by a torus 
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with a large covering factor, and only the reflection from the far side of the torus/disk is 
directly visible. The model assumes the form 

F(E) = [A T T(E, N HtT , T h , E c ) + A R R(E, T h , E c ) + A s E~ Vs + lines]e- N » a » h , (1) 

where T(E, Nh,t, T/i, E c ) is the transmitted emission through an absorption column density 
Nh,t, At is the normalization, is the photon index and Ec the high energy cut-off; 
R(E, Th, E c ) is the reflected nuclear emission by cold dense gas near the central engine or in 
the torus, Ar is its normalization; T s is the photon index of the soft power-law, and A s its 
normalization; lines represents the flux of the emission lines; and Nh is the column density 
of matter external to the nucleus (in the Circinus galaxy and our Galaxy). 

As in § 13.3.11 we consider two cases: the photon-index of the soft power-law is fixed 
to the value from the 0.2-10 keV fit (model 1), and the photon index of the soft power-law 
is tied to the intrinsic power-law index of the nuclear emission (model 2). In model 1 the 
contribution of the soft power-law component to the 13-70 keV PIN count rate is < 5%. 
The column density of matter external to the nucleus Nh is also fixed to the best-fit XIS 
value. The high energy cutoff is fixed to 48.7 keV to avoid unphysical solutions with very 
flat spectra and low cutoff energies. We summarize the best-fit parameters in Table [TJ 

It is obvious that the two models are in fact very similar and the statistics are almost 
identical. Model 2 accounts naturally for the similarity between the soft and hard power- 
law components. Column densities of ~ 4 — 5 x 10 24 are found in both models, which also 
agree with the previous findings from BeppoSax and Integral. In model 1, we estimate the 
luminosity of the nucleus (corrected for absorption) in the 2-10 keV (2-200 keV) band to be 
1.1 x 10 42 ergs -1 ( 2.4 x 10 42 ergs -1 ) assuming a distance of 4 Mpc. 

We finally examine if the spectra can be modeled without the Compton transmitted 
component. In this case, the hard X-ray emission is purely reflected emission. We fit the 
3-70 keV the spectra using the model in § 13.3.11 (case 1). This model yields very flat intrinsic 
power-law with T h ~ 0.34, and unusually large Fe abundance (~ 2.9 solar). We consider 
these parameters unlikely for Circinus galaxy. 

3.3.4- Emission lines and features 

In Fig. we show the ratio of the XIS spectra and our baseline continuum model. The 
best-fit continuum model over 3-70 keV plus the significant emission lines are shown in 
Fig. El Emission lines detected in the 3-10 keV band are listed in Table [2j We detected 
strong features that are consistent emission lines from neutral atoms of Fe and Ni: Fe Ka 
(6.40 keV), Fe K(3 (7.04 keV), and Ni Ka (7.472). A weak line is detected at 6.71 keV, 
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probably from He-like iron Fe Ka. The 7.04 keV line can also be a blend of the "neutral" Fe 
Kf3 and the H-like Fe Ka (6.95 keV and 6.97 keV). We also found a significant line feature 
at ~ 8.2 keV (> 3a). The feature is probably a blend of the H-like Fe K(3 and He-like Fe 
K7 emission lines, another indication of highly ionized gas. However, simply adding a high 
temperature gas component with APEC or MEKAL models in XSPEC cannot improve the 
fit and account for the flux of these lines. All except the 3.74 keV line are unresolved and 
thus modeled with ^-functions to reduce the number of free parameters. 

In Fig. [71 we show the ratio of XIS spectra between 5 and 7 keV and model 1 (§ 13.3.31) 
with only the identified lines shown in Table [2j The spectra show no clear sign of Comp- 
ton Shoulder at 6.23 keV. However, this is likely an artefact due to the inaccuracy of the 
calibration. To illustrate this, we show the XIS spectra of the calibration sources at the 
corners of the XIS detectors (Fig. EJ). The Mn Ka and K/3 lines are modeled with Gaussians. 
The best-fit line width is consistent with 0, but the model line profiles are broader than 
those of the data. We measure the FWHM of the line model convolved with the response 
to be ~ 190 eV, while the FWHM of the line is ~ 180 eV. For the strong Fe Ka line in the 
Circinus galaxy spectra, the broad wings of the response make it hard to properly fit the 
weak Compton shoulder only ~ 0.2 keV below the line energy. In Fig. [HI the broad response 
causes dips ~ 0.1 — 0.2 keV below and above the line energies in the plot of data to model 
ratio. The same is seen in Fig. [51 with a dip at ~ 6.2 — 6.3 keV, which coincident with the 
expected energy of the Compton shoulder. Broad wings are seen in the Fe Ka line, which 
has also been observed in the Chandra observations (Sambruna et al. 2001). Part of the low 
energy wing could attribute to the Compton shoulder. A significant feature is detected at 
~ 5.4 keV. The nature of the line, however, is not known. 

The iron abundance is found to be sub-solar in our models. This seem to disagree 
with the XMM-Newton results (Molendi et al 2003). However, the Fe abundance estimate is 
sensitive to the soft power-law component in our model, which is an over simplification of the 
spectra that includes emission from ULXs, the Thompson scattered nuclear emission, and the 
possible emission from hot gas. If we fit the 3-10 keV XIS spectra without the soft power-law, 
and fix the photon index of the cold reflection component to 1.56, we find Fe abundance to be 
~ 1.35 solar. Therefore, we can not conclude the Fe abundance is actually different from the 
XMM-Newton results. It should be noted that if both reflection and absorption come from 
the torus, the abundance should be the same for the two spectral components. However, our 
absorption model only assumes solar abundance. We estimate how Fe abundance affects the 
absorption column density. In Circinus galaxy, the ion edge is dominated by the reflection 
component, and is determined by the XIS spectrum, while the absorption is determined 
by the PIN spectrum. We fit the PIN spectrum with XSPEC model zvarabs(cutoffpl), 
with the power-law index fixed to 1.56 and the high energy cut-off Ec = 48.7 keV. The 
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abundances of all elements are fixed to the solar value except for ion, which is fixed to the 
best-fit values found above. In this model, if Fe abundance is solar {A Fe = 1), we found 
N H = 266 ± 51 cm" 2 ; If A Fe = 0.75(1.35), the best-fit N H = 318 ± 60(216 ± 41) cm" 2 . In 
either case, the under or over estimate of N H by assuming A Fe — 1 is at ~ la level. 

If the flux in the 7.04 keV line is from Fe K/?, the flux ratio of Fe K(3 and Ka would be 
0.20 ±0.02. This value has not been corrected for the possible loss due to the low energy tails 
of the lines and the Compton Shoulder. The result is consistent with that from Chandra 
(0.22 ± 0.05, Sambruna et al. (2001) and BeppoSax (0.191, Guainazzi et al. (1999), but 
appear to be significantly higher than that from XMM-Newton (0.14 ± 0.01, Molendi et al. 
2003) and the expected value ~ 0.16 (Basko 1978). If the flux in the 7.04 keV is a blend of 
the "neutral" Fe Kf3 and the H-like Fe Ka, which we think likely based on our detection of 
lines that are consistent with the highly ionized species, the different line ratio between the 
XMM-Newton and other observations may suggest variability in the ionization state of the 
gas, and could potentially constrain the physical size of the ionization region. The detailed 
study of the ionized gas in Circinus galaxy is beyond the scope of this paper. 

4. Concluding remarks 

Our main results from an X-ray observation of the Circinus galaxy with the Suzaku 
satellite are: 

1. There is minimal contamination (±5%) of the 0.2-10 keV radiation by the periodic 
ULX CG X-l, and less at higher energies. 

2. Below 10 keV, we see a soft power-law ( that could come from a combination of extended 
emission, contaminating point sources in the galaxy, or the Thompson scattered nuclear 
emission) and a strong cold reflection component (as observed by earlier satellites) from 
the nuclear emission. 

3. The spectrum above 20 keV is well described by a heavily absorbed Compton thick 
(N H = (4 — 5) x 10 24 cm" 2 ) power- law model with a high energy cut-off (E c ^ 49 keV). 
Based on the lack of variability in the PIN light curve, we suggest that the nucleus 
is obscured by a torus with a large covering fraction, with the reflection component 
originating from the far side of the torus. 

4. The galaxy is detected in the 50-200 keV band at 2.5a level. This flux agrees well 
with an extrapolation to the higher energies of the transmitted component seen in the 
13-50 keV band. 
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5. We have detected strong emission lines at energies consistent with the florescent lines 
of neutral Fe Ka, K/3, and Ni Ka. These features are consistent with the reflection 
from cold dense gas in the torus. However, we were not able to find the Compton 
shoulder in our XIS spectra. We point out this non-detection can be attribute to a 
calibration problem. 

6. We also found emission features which are consistent with emission lines from He- 
like Fe Ka (6.71 keV), H-like Fe K/3 (8.2 keV). The emission line at 7.04 keV is also 
consistent with the line energy of H-like Fe Ka. These suggest that highly ionized gas 
is present in the Circinus galaxy. 

Y. Yang would like to thank Dr Richard Mushotzky for discussions on the XIS spectra. 
We thank Drs. Koji Mukai, Diana Harms, and Kenji Hamaguchi at Suzaku GOF for help 
with issues related to data reduction and calibration. We thank the anonymous referee for 
suggestions that improve this paper. This work is supported by NASA through LTSA grant 
NAG 513065 and Suzaku NMX06AI34G granted to the University of Maryland. 
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Table 1. Best-fit Continuum 





Model 1 


Model 2 


N h .t (10 22 cm- 2 ) 




470l^° 




u - lo -0.04 


14+ 004 


r h 


i cc+0.09 


1 ro+0.07 

1.08_ 10 


E c (keV) 


48.7 (fixed) 


48.7 (fixed) 


Ar & 


(2.2l°; 2 ) x 10- 2 


(2.4l°' 2 ) x 10- 2 


Ap e (solar) 


0.75 ±0.09 


74+ - 10 

u - '^-0.06 


T s 


1.54 (fixed) 


1.58 = T h 


A s a 


1.1 x 10- 3 (fixed) 


(1.0 ±0.3) x 10- 3 


N H (10 22 cm- 2 ) 


0.38 (fixed) 


0.38 (fixed) 


X 2 /d.o.f h 


2066/1814 


2067/1813 



a unit: phkeV 1 cm 2 s 1 at 1 keV. 

b The x 2 are obtained with added emission lines listed in Ta- 
bled 



Table 2. Emission lines 



Energy (keV) Identification flux a Gaussian a (keV) b EW (keV) 



3.12 ±0.02 


Ar XVII 


7.2 + f° 7 x 10- 6 





0.032 


3.75 ±0.05 


Ca XIX, Ar XVII 


7.8ll°';°| x 10- 6 


(7.3 ±0.05) x 10- 2 


0.040 


5.38 ±0.04 


unkown 


4.22 x 10- 6 





0.022 


6.405 ±0.001 


Fe Ka (II-XVII) 


(3.45 ±0.04) x 10- 4 





1.49 


6.716 ±0.006 


Fe Ka (>XVII) 


(4.7 ±0.02) x 10- 5 





0.049 


7.040 ± 0.004 


Fe K/3 


(6.95 ±0.02) x 10- 5 





0.126 


7.48 ±0.01 


Ni Ka 


(1.66 ±0.02) x 10- 5 





0.081 


8.28 ±0.04 


Fe K/3 (XXVI)+ K 7 (XXV) 


6.01 x 10- 6 





0.058 



a Unit: ph cm 2 s 1 

b Most of the lines have a consistent with 0, so their width is set to to reduce the number of free 
parameters. 
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Fig. I. — Chandra 0.5-7 keV ACIS-S image of the Circinus galaxy (Obs Id 356, PI: Andrew- 
Wilson). The image is slightly smoothed for better viewing. The 2.5' radius circular XIS 
spectral extraction region is shown. The two brightest ULXs included in the extraction 
region, CG X-l, and X-2 in Bauer et al. (2001), are marked. A bright point source ~ 5' 
SW of the Circinus galaxy is also seen in the Chandra image. This source is within the 
field-of-view of the HXD detector. 
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Fig. 2. — The folded 0.2-10 keV XIS light curve of the Circinus Galaxy. Light curves from 
all of the four XIS instruments are included. The folding period is 27 ks and the count rate 
is normalized to the mean value of 2.25 cts/s. 
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Fig. 3. — The 3-200 keV Suzaku spectrum of Circinus galaxy and the folded best-fit model 
(model 1 in the text). The color representation for the spectra from each of the instruments 
is defined as: black-XISO, red-XIS2, green-XIS3, blue-HXD/PIN, and cyan-HXD/GSO. 
Note that while the GSO spectrum agree very well with the extrapolation of the best-fit 
XIS+PIN model, it is not included in the spectral models. 
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Fig. 4. — x 2 as a function of I\ and Eq when the 20-70 keV PIN spectrum is modeled 
with a nuclear power-law viewed through a Compton transmitting medium (see text). The 
contours correspond to 1, 2 and 3-sigma confidence levels. 
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Fig. 5. — The ratio of the XIS spectrum vs. the best-fit continuum. The spectra have been 
rebinned for the purpose of viewing only. Color representation is the same as in Fig. [3J 
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Fig. 6. — The unfolded spectrum between 3-70 keV band plotted in vF v . Model 1 in the 
text has been used in unfolding the spectra. Color representations are the same as those in 
Fig.S 
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Fig. 7. — The ratio of the XIS spectrum and the model with the identified lines. The spectra 
have been rebinned for the purpose of viewing only. Color representation is the same as in 
Fig.S 
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Fig. 8. — Upper panel: The spectra of the calibration sources in XIS-0,2,3. The best-fit 
5-function models folded with the responses are also shown. Lower panel: the ratio of the 
data and models in the upper panel. The color scheme is the same as Fig. [3] 



